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The DGP brane-world model provides a simple alternative to the standard LCDM cosmology, 
with the same number of parameters. There is no dark energy - the late universe self-accelerates 
due to an infrared modification of gravity. We compute the joint constraints on the DGP model 
from supernovae, the cosmic microwave background shift parameter, and the baryon oscillation peak 
in the SDSS luminous red galaxy sample. Flat DGP models are within the joint 2 sigma contour, 
but the LCDM model provides a significantly better fit to the data. These tests are based on the 
background dynamics of the DGP model, and we comment on further tests that involve structure 
formation. 

PACS numbers: 95.36.+X, 98.80.-k, 98.80.Es 



1. THE DARK ENERGY PROBLEM 

The acceleration of the late-time universe, as implied 
by observations of supernovae, cosmic microwave back- 
ground anisotropics and the large-scale structure, poses 
one of the deepest theoretical problems facing cosmol- 
ogy Q- Within the framework of general relativity, the 
acceleration originates from a dark energy field (or ef- 
fective dark energy) with negative pressure {w = p/p < 
— |), such as vacuum energy {w = —1) or a slow-rolling 
scalar field ("quintessence", w > —1). So far, none of the 
available models has a natural explanation. 

For the simplest option of vacuum energy, i.e., the 
LCDM model, the incredibly small value of the cosmo- 
logical constant 
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cannot be explained by current particle physics. In ad- 
dition, the value needs to be fine-tuned, 
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which also has no natural explanation. Quintessence 
models attempt to address the fine-tuning problem, but 
do not produce a natural solution - and also cannot ad- 
dress the problem of how A is set exactly to 0. 

Alternatively, it is possible that there is no dark en- 
ergy, but instead an infrared modification of general rel- 
ativity, i.e., on very large scales, r > Hq^, that accounts 
for late-time acceleration. (Note that this does not re- 
move the problem of explaining why the vacuum energy 
does not gravitate.) Schematically, we are modifying the 
geometric side of the field equations. 
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rather than the matter side, 

Gfj,u — SttG {Tf^i, + T^^''^) 



as in general relativity. It is important that the modifica- 
tion is covariant and incorporates deviations from homo- 
geneity and isotropy, so that one can compute not only 
the background dynamics, but also the perturbations. 



2. DGP MODIFIED GRAVITY 

One of the simplest covariant modified-gravity models 
is based on the Dvali-Gabadadze-Porrati (DGP) brane- 
world model, as generalized to cosmology by Deffayet 
(It is worth noting that the original DGP model with a 
Minkowski brane was not introduced to explain acceler- 
ation - the generalization to a Friedman brane was sub- 
sequently found to be self-accelerating.) In this model, 
gravity leaks off the 4-dimensional brane universe into 
the 5-dimensional bulk spacetime at large scales. At 
small scales, gravity is effectively bound to the brane 
and 4D gravity is recovered to a good approximation, 
via the lightest modes of the 5D graviton - effectively 
via an ultralight metastable graviton in the 4D universe. 
The transition from 4D to 5D behaviour is governed by a 
crossover scale r,,', the weak- field gravitational potential 
behaves as r~^ for r <^ Tc and as for r ^ rc- 

The energy conservation equation remains the same as 
in general relativity, but the Friedman equation is mod- 
ified: 
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These equations imply (for the CDM case p = 0) 
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Equation ((T)) shows that at early times, when + 
K/a^ ^ r~^, the general relativistic Friedman equation 
(5) is recovered. By contrast, at late times in a CDM uni- 
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verse, with p oc a ^ 0, we have 

H Hnc, — - 



(9) 



Gravity leakage at late times initiates acceleration - not 
due to any negative pressure field, but due to the weak- 
ening of gravity on the brane. Since Hq > i/oo, in order 
to achieve self-acceleration at late times, we require 



(10) 



and this is confirmed by fitting observations, as discussed 
below. 

In dimensionless form, the modified Friedman equa- 
tion ^ is 
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From Eq. ||SJ, the dimensionless acceleration is 

2VLr^ - r2m(l + 2)^ 
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so that the redshift when acceleration starts is given by 
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The (i7,.^,f2iu) plane is illustrated in Fig.Q] This shows 
the flat model curve [from Eq. H12|l ] , 
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and the present-time zero-acceleration curve [from 
Eq. CHI)], 
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By checking for real roots z > of the equation H{z) = 0, 
given by Eq. (|ll|l with radiation included, we also find 
the region of the plane containing bouncing models, i.e., 
models that are currently expanding but were collapsing 
in the past. 

The modified Friedman equation in DGP may reinter- 
preted from a standard viewpoint. We define the effective 
dark energy density poff = SH/SnGrc- Then the effective 
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FIG. 1: The parameter space for DGP models. The dotted 
(red) curve is the flat case, with closed models above and open 
models below; the dot-dashed (black) curve demarcates mod- 
els that are currently accelerating (above) from those decel- 
erating (below); the shaded (blue) region contains expanding 
models that bounce in the past, i.e. do not have a big bang; 
the likelihood contours are fits to the SNe Gold data [dashed 
(brown)] and Legacy data [solid (blue)]. 



dark energy equation of state WeS = Pos/PcS is given by 
Poff + 3if(l -I- Weff)peff = 0. Thus Peff and Wcft give a 
standard general relativistic interpretation of DGP ex- 
pansion history, i.e., they describe the equivalent general 
relativity dark energy model. For the flat case, Hk = 0, 
we find 
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which implies 
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The DGP and LCDM models have the same number of 
parameters, with substituting for A. (Neither model 
provides a natural solution to the late-acceleration prob- 
lem - both A and Tc must be fine-tuned to match ob- 
servations.) The DGP cosmology therefore gives a very 
useful framework for comparing the LCDM general rela- 
tivistic cosmology to a modified gravity alternative. Does 
the DGP model pass the observational tests? The next 
section considers this question. 



3. OBSERVATIONAL CONSTRAINTS ON DGP 

The fundamental test of the background dynamics of 
a cosmological model is the SNe magnitude- redshift test, 
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FIG. 2: Likelihood contours for the CMB shift parameter S 
in DGP (left) and LCDM (right). 



based on the luminosity distance, 
(1 + ^) 
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where 



T{x) = (a;,sina;,sinh2;) for is: = (0, 1,-1). (21) 

The 68, 95 and 99% likelihood contours from fits to the 
SNe data 0, 3 are shown in the DGP parameter plane 
in Fig. HI 

The luminosity distance is determined purely by the 
modified Friedman equation (|lll) . and there is a general 
relativistic dark energy model that can exactly mimic 
the DGP luminosity distance. This equivalent model has 
dark energy equation of state equal to the DGP effec- 
tive equation of state Wcs, as given in Eq. lfTH|) for the 
flat case. Of course, there is no physical justification 
for the equivalent general relativistic model, whereas the 
DGP model has a clear physical motivation. Neverthe- 
less, tests based on the Hubble rate H{z) cannot dis- 
tinguish the DGP from a GR dark energy model. Fur- 
ther tests are needed to discriminate DGP from general 
relativistic models. The CMB anisotropics and matter 
power spectrum provide in principle suitable discrimina- 
tory tests. These tests require a detailed understand- 
ing of the evolution of density perturbations in the DGP 
model. Unfortunately the problem is highly complicated 
and not yet solved, essentially since density perturba- 
tions are coupled to perturbations in the 5D gravitational 
field |5|. 

Previous results on the CMB anisotropics and struc- 
ture formation in DGP provide valuable strategies for 
testing DGP and discriminating it from general relativis- 
tic models. However, these results are based on neglect- 
ing the 5D perturbations, which effectively imposes an 
inconsistency, leading to a violation of the 4D Bianchi 
identity, as explained by Koyama and Maartens Q. On 
sub-Hubble scales, linear density perturbations in the 
DGP model have been computed in a way that satisfies 
the 5D equations Q , and this result confirms and gener- 
alizes the results of Lue and Starkman (3] ■ Further work 



is needed to quantify the corrections to general relativis- 
tic perturbations on scales near and above the Hubble 
radius. 

Until the problem of cosmological perturbations in 
DGP is solved, we must turn to other observations, inde- 
pendent of the SNe redshifts, that rely on the background 
model and do not rely (or rely only weakly) on the den- 
sity perturbations. 



CMB shift 



The CMB shift parameter 
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relates the angular scale of the first acoustic peak to the 
angular diameter distance to last scattering and the phys- 
ical scale of the sound horizon. It is effectively model- 
independent and insensitive to perturbations, and may 
be used to constrain the background dynamics of models, 
such as quintessence models Q. Wang and Mukherjee 
have used the WMAP 3-year data to derive 



5 = 1.70 ±0.03, 
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with Zr = 1090. The likelihood contours for S in the 
DGP and LCDM cases are shown in Fig. [3 Note that 
the DGP contours are more in the open-model region for 
fim < 0.4 than the LCDM contours. 



Baryon oscillations 

Fairbairn and Goobar used the baryon acoustic 
oscillation peak recently detected in the SDSS luminous 
red galaxies (LRGs) jlj as an independent test of the 
DGP. The correlation function for SDSS LRGs shows a 
peak at a scale ~ 100/i~^ Mpc, corresponding to the first 
acoustic peak at recombination (determined by the sound 
horizon). The observed scale effectively constrains the 
quantity [l^ 
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and Eisenstein et al. find that ll^ 



A = 0.469 ±0.017, 



(24) 



(25) 



where zi = 0.35 is the typical LRG redshift. (Following 
Ref. |ll| , we have suppressed a weak dependence of A on 
the spectral tilt.) 

We should point out that there is a level of uncer- 
tainty in the use of the BO measure A. The BO data 
are analyzed via a fiducial LCDM model, with a single 
scaling relation used to make a best fit, and a compres- 
sion of the data to a constraint at a single redshift [T^ . 
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FIG. 3: Joint constraints [solid thick (blue)] on DGP models from the SNe data [solid thin (yellow)], the BO measure A [dotted 
(green)] and the CMB shift parameter S [dot-dashed (red)]. The left plot uses SNe Gold data, the right plot uses SNLS data. 
The thick dashed (black) line represents the flat models, fiic = 0. 




FIG. 4: As in Fig. but for LCDM. 



Within the LCDM class of models, this approximation is 
accurate to within a few percent 12] (less than current 
errors in the determination of Jim). As pointed out by 
Dick et al. this single-redshift feature, which may 
be robust for LCDM and constant-w models, could in- 
troduce significant errors for models with nonconstant 
w. The DGP has an effective nonconstant w, although 
over the SDSS LRG redshift range (0.16 < z < 0.47) 



the variation is small. It is not clear whether the errors 
introduced in fitting DGP models are at the acceptable 
level of a few percent. It is effectively implicitly assumed 
that this is the case in Refs. |Tll ITl . Hsj l. This may be 
a reasonable assumption, since the DGP comoving lumi- 
nosity distance is close to LCDM in the redshift range 
for the LRGs 0. Ultimately, the BO data should be re- 
analyzed without assuming a fiducial LCDM model and 
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without compression to a single redshift, so as to be a 
more reliable constraint on alternatives to LCDM. 

Another issue is the extent to which the BO con- 
straint relies on the LCDM matter power spectrum. The 
DGP corrections to the shape of the matter power spec- 
trum at last scattering on scales around the first acoustic 
peak scale are small, although there may be some scale- 
dependence in the corrections 16\. 

Fairbairn and Goobar showed that the joint constraints 
from SNe (Legacy) data and the BO measure A ruled out 
a flat DGI' model at 3a. Alam and Sahni used the 
Gold SNe data with the BO measure A and showed that 
the flat DGP model is marginally allowed at 2a. (See 
also Refs. A conservative conclusion from these 

results is that the SNe-BO constraints exclude flat DGP 
models at 2a, unlike the case of LCDM models, where 
flat models are within the joint la contour 

However, this negative conclusion is over-turned when 
we include the CMB constraint from S. 



Combined SNe-CMB-BO constraints 

The effect of the CMB shift S is to puU the best-fit 
model down from the closed region towards the flat curve. 
This is clearly illustrated in Fig.|21 without the shift con- 
straint, the joint SNe-BO best-fit would be well above the 
flat model curve, as found by Fairbairn and Goobar [T]| . 

The corresponding constraints for the LCDM model 
are shown in Fig. 21 and are consistent with the results 
of Ichikawa and Takahashi 17] . For LCDM, the SNe- 
BO constraints already allow a flat model at la, and 
the effect of the shift parameter is simply to narrow and 
rotate the contours, ft is clear from Figs. |21 and ^ that 
there is some tension between the data and DGP, which 
is not the case for LCDM. 



is signalled by the rather large value of Wcs{0): using 
f^,„ « 0.26 from Tables I and II, Eq. lO gives 
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best-fit 
acceleration 
parameter 


best-fit 
density 
parameter 


best-fit 
curvature 
parameter 


value 


DGP 
LCDM 


flr^ = 0.125 

= 0.730 


flm = 0.270 
= 0.285 


fix = +0.0278 
flK = -0.0150 


185.0 
177.8 


TABLE I: Best-fit parameters from SNe (Gold)-CMB shift- 
BO constraints, and values, for the DGP and LCDM mod- 
els. 




best-fit 
acceleration 
parameter 


best-fit 
density 
parameter 


best-fit 
curvature 
parameter 


value 


DGP 
LCDM 


ilr, = 0.130 

= 0.740 


fin, =0.255 
fin, = 0.270 


^A' = -h0.0300 
^A' = -0.0100 


128.8 
113.6 



TABLE II: As in Table I, for the Legacy SNe data. 



For the Gold data, LCDM reduces the by 7, 
whereas the Legacy data gives an even larger reduction 
of 15. LCDM provides a significantly better fit to SNe- 
CMB shift-BO observations than DGP, especially for the 
Legacy SNe data. 

The P6r degree of freedom provides a comparison 
of the goodness of fit across the two SNe data sets: 



4. CONCLUSIONS 

The combination of three independent observational 
constraints - SNe, CMB shift parameter, baryon oscil- 
lations - is a powerful way to test the background dy- 
namics of a cosmological model. It has been used for 
quintessence models by Ichikawa and Takahashi ^3] and 
for an f{R) modified gravity model by Amarzguioui et 
al. We applied this triple joint constraint to the 

DGP case, using both the Gold and the Legacy data, 
and showed that flat DGP models are within the 2a joint 
contour. We note that the best-fit DGP model is slightly 
open, as opposed to the slightly closed nature of the best- 
fit LCDM model. 

Although flat DGP and LCDM models are both within 
the 2a contour, it is clear from Figs. El and 01 that LCDM 
fits the data more easily. A quantitative measure of this is 
via the values for the best-fit models in the two cases, 
as shown in Tables I (Gold SNe data) and II (Legacy 
SNe data). The tension between the data and the DGP 



(xVdof)j,Qp 



(xVdof) 



LCDM 



1.179 Gold 

1.120 Legacy 

1.132 Gold 

0.988 Legacy 



The Legacy data leads to a better fit for both DGP and 
LCDM. 

One qualification to our results is that we have not 
shown the reliability of the BO measure A as a clean test 
of the background dynamics for DGP (this also applies 
to previous applications of the BO constraint to DGP 
models 

[H III El)- The BO measure A is tailored to 
the LCDM model, as explained above. Although we may 
expect that only small corrections are involved, it is not 
clear whether this introduces an artificial bias against the 
DGP model. This is a nontrivial problem to resolve, and 
further work is needed. 

At a more fundamental level, the DGP needs to be 
tested against the data on CMB anisotropics and the 
matter power spectrum, and these tests provide also the 
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means to discriminate DGP from the equivalent general 
relativistic dark energy model. However, as explained 
above, theses tests cannot yet be carried out, since the 
analysis of the density perturbations in DGP has not yet 
been sufhciently developed @| . We should also note that 
there is a ghost at the linear perturbative level in the 
asymptotic de Sitter state of the self-accelerating DGP 
cosmology jl^. Classically, this ghost is not a problem, 
but it does render the quantum vacuum unstable, and we 
effectively assume that an ultraviolet completion of the 
DGP model may be found that cures the ghost instabihty 
in the quantum vacuum. 
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